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O ' Abstract: 

O . 

CN . The critical behaviour of many spin models can be equivalently formulated as perco- 

lation of specific site-bond clusters. In the presence of an external magnetic field, such 
clusters remain well-defined and lead to a percolation transition, even though the system 
no longer shows thermal critical behaviour. We investigate the 2-dimensional Ising model 
and the 3-dimensional 0(2) model by means of Monte Carlo simulations. We find for 
^ , small fields that the line of percolation critical points has the same functional form as the 

>• [ line of thermal pseudocritical points. 

o 

CO 

Percolation theory [[l], [2[ provides in many cases an elegant interpretation of the mech- 
anism of second order phase transitions: the existence of an infinite spanning cluster 
represents the new order of the microscopic constituents of the system due to sponta- 
neous symmetry breaking. 

In the Ising model, a rigorous correspondence between thermal critical behaviour 
and percolation was established by Coniglio and Klein ||; the clusters are constructed 
by linking like-sign spins through a temperature-dependent bond probability p = 1 — 
exp(—2J/kT), where J is the Ising coupling and T the temperature. This result has been 
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recently extended to a wide class of models, from continuous spin Ising-like models 



to 0(n) models ||. In all cases the system was considered in the absence of any external 
field. The reason for this is clear: by introducing an external field H, we explicitly break 
the symmetry of the Hamiltonian of the system, thus eliminating the thermal critical 
behaviour of the model. None of the thermodynamic potentials exhibits discontinuities 
of any kind, since the partition function is analytical for H ^ 0. 

On the other hand, the Coniglio-Klein clusters can be built as well when H^O. Because 
of the field, the system has a non-vanishing magnetization m parallel to the direction of 
H for any finite value of the temperature T. For T — * oo, m— > 0, and for T = 0, m = 1. 
This suggests that for a fixed value of H, the clusters will start to form an infinite network 
at some temperature T P (H). Varying the field H, one thus obtains a curve T P (H) in the 
T — H plane, the Kertesz line JF]. 

The Kertesz line specifies the usual percolation threshold, just as in the case H = 0, 
and the percolation variables exibit the usual singularities, leading to a set of critical 
exponents. In particular, the percolation strength, which is the relative size of the perco- 
lation cluster compared to the total volume of the system, remains the order parameter 
of the percolation transition. 



The existence of such a line of genuine critical behaviour also for H^O could provide 
a criterium to define different phases and the transition between them in a more general 
sense ||. For this purpose, it is important to find out whether there is any relationship 
between the geometrical singularities at the Kertesz line and thermal properties of the 
system along this line; apparently, rather little is known about such connections |9fl-fTT| . 



In the present work we want to address this aspect in more detail. In particular, we 
shall investigate the relation between the Kertesz line and the so-called "pseudocritical" 
or " crossover" line, defined by the values of the temperature at which the magnetic sus- 
ceptibility peaks as function of the field. In general, the two lines cannot coincide, since 
when H — > oo, the Kertesz line leads to a finite T, while for the pseudocritical line this 
results in T — > oo; we return to this point later on. We shall here study two spin models, 
the Ising model in two dimensions and the 0(2) model in three dimensions, and calculate 
in both cases the functional form of the Kertesz line and the pseudocritical line in the 
limit of a small external field. 

The Hamiltonian of the models at study can be written in the following general form: 

J H 
= -p^E s i s j-^E s HSi, (1) 

(ij) i 

where Si, Sj are two-dimensional unit vectors for 0(2) and simple scalars (±1) for the Ising 
model; similarly, sh is a unit vector (or a scalar) in the direction of the external field. We 
now first consider the percolation problem. 

In Fig. [1] we show schematically the Kertesz line for the Ising model. For a vanishing 
external field, H = 0, the Coniglio-Klein cluster definition with 

PlJ = p b = 1 - exp{-2J/A;T} (2) 

as bond weight between two adjacent spins i,j assures that we recover the usual thermal 
threshold of the Ising model. For H — > oo, all lattice spins will be aligned with the 
field at any temperature T. However, the bonds between adjacent spins will be occupied 
according to the bond weight p b , so that now the site-bond problem turns into pure 
bond percolation. The percolation transition will then take place for that value of the 
temperature T b for which the probability pij equals the critical density pb{d) of random 
bond percolation in d dimensions, 

p b (d) = 1 - exp{-2J/A;T b }, (3) 

i.e., for 

(4) 



k log[l -Pb(d)\ 

These arguments can simply be repeated for the Kertesz lines corresponding to 0(n) spin 
models. 

As mentioned above, each point of the Kertesz line is a standard percolation point. 
At H = 0, the critical exponents of the percolation variables coincide with the thermal 
critical exponents of the magnetization transition (Ising, 0(ri)). However, for any H^O, 
they are predicted to switch into the (different) universality class of random percolation 
in the same dimension. We will verify this for each of the percolation points we shall 
determine. 
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Figure 1: Schematic behaviour of the Kertesz line for the Ising model 



Turning now to the pseudocritical thermal behaviour of the models considered, we 
define the reduced temperature t and external field h measured relative to the spin-spin 
coupling strength 



t 



T — T r . 



H 



Tc • ■- J < 5 > 

For spin systems, the renormalization group approach leads to the following expression 
of the magnetization m when t, h <C 1 



m = h^fcit/h 1 ^), 



(6) 



where fa is a scaling function and /3, S are the critical exponents for the magnetization 
transition at h — 0. We can rewrite Eq. (M) in the form 



(7) 



with a new scaling function Fq] the entire dependence on the field h is now put into the 
function F G . 

The susceptibility x is given by 



X - (0) oc t^F' G {x)th-^ s -\ 



(8) 



where F G (x) is the derivative of the function Fq with respect to its argument x = t/h 1 ^ 6 . 
This time we move the dependence on the reduced temperature t into a new function T 
of the argument x and obtain 



X 



(9) 



To find the equation of the pseudocritical line, we have to determine the temperature t x (h) 
at which the susceptibility x peaks for a given value h of the field. Hence the derivative 
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of x with respect to the reduced temperature t must vanish at t = t x , 

= o - hV'-WWl^a/hW*) = /i 1/5 -^ 1//35 r(x)| t=tx = 0; (10) 



dx 
dt 



here F'(x) \t=t x is the derivative of the function T with respect to its argument x calculated 
at t = t x . The derivative dx/dt\ t =t x can be zero only if !F'{x)\t=t x = 0. That will occur 
at some value x x of the argument x = t x //i 1///3<s , which provides the relation between t x 

tjh 1/ps = x x — t x oc h l/l3& . (11) 

The procedure is obviously independent of the value of the field h, so that as long as h is 
small, the pseudocritical line is described by a simple power law. 

We stress that, for h — > oo, the temperature of the susceptibility peak diverges, t x — > 
oo, whereas we have seen that the Kertesz line has a finite endpoint, given by Eq. @. 
Hence the two curves, which start from the same point but tend to two different limits, 
must certainly differ for sufficiently large values of h. We want to investigate here, however, 
what happens in the vicinity of the thermal critical point of the model, i.e. for small fields 
(hmO). 

In order to study the percolation transition we need to redefine each given spin con- 
figuration as a cluster configuration, by grouping all spins of the lattice in clusters. We 



do this with the help of the algorithm of Hoshen and Kopelman |fL2|l , using free bound- 
ary conditions for the cluster identification and assuming that a cluster percolates if it 
connects each pair of opposite sides (faces in 3d) of the lattice. 

After measuring the size of all clusters, i.e., the number of sites in each cluster, we 
evaluate the following two quantities: 

• The average cluster size S, 

S = , (12) 

Es n s s 

where n s is the number of clusters of size s; the sums exclude the percolating cluster. 

• The percolation strength P, 

p size of the percolating cluster 
total number of lattice sites ' 

which is the order parameter of the percolation transition. 

In the infinite volume limit, the percolation variables are described by simple power laws 

P ot (T p -Tfp (T<T P ), S<x\T-T p \t (14) 

sufficiently close to the critical threshold T p . 

In order to determine the critical point of the percolation transition, we exploit the 
properties of a further very useful variable. For finite lattice size, there can be spanning 
clusters even for temperatures T > T p , and there are spin configurations at temperatures 
T < T p below the critical threshold without such clusters. The probability of finding a 
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spanning cluster on a finite lattice of linear dimension L at a temperature T is a well- 



defined function, 
behaves as 



the percolation cumulant U(T,L) [13[]. For T ~ T p and large L, it 



I1(T, L) = $ 



T-T P 



(15) 



where z/ p is the critical exponent governing the divergence of the percolation correlation 
length. The function I1(T, L) is not a genuine percolation variable, since it has a non- 
trivial meaning only on finite lattices. On an infinite lattice it reduces to a step function: 
it is zero for T > T p and unity for T < T p . Nevertheless, its particular features make 
it a powerful tool to extract information about the critical properties of percolation. In 
particular, at T = T p , H(T p , L) = $(0) for any value of L. That means that if we calculate 
the percolation cumulant as a function of T for different lattice sizes, all curves IT(T, L) 
will cross at the critical temperature T p . Moreover, if we plot the different IT(T, L) as 
a function of the variable X = [(T — T p )/T p ]L 1 ' Vp , the resulting functions $(Jf) must 
coincide for all lattice sizes. To define the scaling function <&(X), one also has to know 
the value of the exponent v p . We can therefore determine v p by looking for the best 
scaling form of the percolation cumulant curves. 

Let us now first consider the results obtained for the 2D Ising model. Our aim is to 
determine the functional form t p (h) describing the Kertesz line. In order to compare it 
to the pseudocritical line, we have to perform simulations at small values of the reduced 
field h. We consider the following four values: h= 0.002, 0.001, 0.0005 and 0.00025. The 
Monte Carlo algorithm used in the simulations is the Wolff cluster update extended to 
the case of a non-vanishing external field 



HI. For each iteration we have measured the 



energy e of the configuration, its magnetization m, and the percolation variables S and 
P. To determine the percolation critical point with a good accuracy, we work with rather 
large lattices: 300 2 , 500 2 , 700 2 , 1000 2 . This is helpful also for another reason. Since the 
field is quite weak, the scaling behaviour of both thermal and percolation variables could 
be perturbed by the "near-by" case of the system without field, which could affect the 
values of the critical exponents. Such perturbation becomes smaller the larger the size of 
the lattice. 
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Figure 2: Percolation cumulants for four lattice sizes; (3 = J/kT, h=0. 00025. 



Fig. shows the percolation cumulants as a funtion of (3 = J/ kT for different lattice 
sizes, at the smallest field we have studied, h = 0.00025. The curves cross clearly at the 
same point (fi p = 0.43933). To check the universality class of the percolation exponents in 
this case, we rescale the percolation cumulants as discussed in Section 3. Using J/kT p = 
0.43933, we consider two options for the critical exponent u p , the 2-dimensional random 
percolation value vpp = 4/3 and the 2D Ising exponent ujm = 1. The results are shown 
in Figs. |^ and they clearly indicate that the critical behaviour at the Kertesz line is 
determined by the random percolation exponent vrp, as expected. 
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Figure 3: The rescaled percolation cumulants of Fig. [ 
2-dimensional random percolation exponent vpp = 4/3. 
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Figure 4: The rescaled percolation cumulants of Fig. || with J/kT p = 0.43933 and the 
2-dimensional Ising exponent ujm = 1- 
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Using a standard finite-size scaling analysis of the average cluster size S near the 
critical point, we obtain for the ratio of the exponents 7 p /f p = 1.80(2), in accord with the 
value predicted from the random percolation exponents (129/72=1.7916) and not with 
that from the Ising model (7/4=1.75). Repeating the percolation cumulant study for the 
other three values of h, we find the critical percolation temperatures listed in Table [I] and 
plotted in Fig. |5[ In the Figure, the reduced percolation temperature t p = (T p — T c )/T c 
is shown as function of the external field; here T c is the critical temperature of the Ising 
model without field (J/kT c = (1/2) log(V2 + 1) = 0.44068679). On the logarithmic scale 
of the plot the data points fall remarkably well on a straight line. We thus conclude that, 
for the Kertesz line of the 2D Ising model, 



t p = a p h K , for/i— >0, (16) 

where a p is a constant of proportionality. From the slope in Fig. [| we obtain k = 0.534(3). 
On the other hand, the exponent of the thermal pseudocritical line is 1//35 (see Eq. fllTD ). 
For the 2D Ising model, j3 = 1/8, 5 = 15, so that 1/(35 = 8/15 = 0.533. The agreement 
between this value and the exponent found for the Kertesz line is excellent. 



h 


J/kT p 


0.00025 


0.43933(3) 


0.00050 


0.43875(4) 


0.00100 


0.43786(4) 


0.00200 


0.43662(4) 



Table 1: Percolation thresholds for the 2D Ising model at different values of the field. 
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Figure 5: Kertesz line and pseudocritical line for the 2D Ising model. 




Figure 6: Susceptibility x for the 700 2 and 1000 2 lattices at h = 0.001. 



We thus ask ourselves whether the Kertesz line and the pseudocritical line overlap 
in the range of h values we have chosen. To answer this, we determine some points of 
the pseudocritical line and compare them to the corresponding points of the Kertesz line. 
Fig. | shows the susceptibility as a function of (3 for h = 0.001, calculated on the two 
largest lattices we have studied. The curves are obtained by interpolating the data points 



using the density of state method The Figure clearly indicates that there is no 

noticeable change of the peak between 700 2 to 1000 2 , so that we are effectively at the 
infinite volume limit. The position of the peak in the Figure, however, clearly differs from 
the corresponding percolation point at h = 0.001, indicated by the dotted line. 

To complete the comparison, we determine the positions of the susceptibility peaks 
for two other values of the field, h = 0.0005 and h = 0.002. As above, we find clear 
discrepancies with the corresponding points of the Kertesz line; the results are included 
in Fig. |5|, where it is also seen that the power law of Eq. ( |TT| ) 

t x = a x /i 1///3<5 , a x = const. (17) 

continues to be valid with great precision. We conclude that for the 2D Ising model the 
functional dependence of the reduced temperature t on the field h for h— >0 is for the 
Kertesz line and for the pseudocritical line described by the same function, determined 
by the thermal critical exponent 1/(38; but the two curves do not coincide, since a p ^ a x , 
as is evident in Fig. |^. 

We now turn to the corresponding study of the 3D 0(2) model. Here we have deter- 
mined four points of the Kertesz line, for external field h values 0.001, 0.00050, 0.00025, 
and 0.000125, using again the cluster update |L4] as Monte Carlo algorithm, on lattices 



of size 30 , 40 , 50 , 60 , 70 , 80 , 100 . The clusters we considered are those used in the 



Wolff algorithm with no external field (see || |T6||). 

The positions of the percolation points, determined as before by the crossing of the 
percolation cumulant curves, are listed in Table |2j. As for the Ising model, we have checked 
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that the universality class of the exponents of each percolation transition is that of random 
percolation, here in three dimensions. We show these points in Fig. |7|, again plotting the 
/i-dependence of the reduced percolation temperature t p = (T p — T c )/T c . The critical 
temperature T c of the magnetization transition for the 3D 0(2) model without field was 
found to be J/kT c = 0.454165 in |7|. 



h 


J/kT p 


0.000125 


0.45312(2) 


0.000250 


0.45257(2) 


0.000500 


0.45178(3) 


0.001000 


0.45053(3) 



Table 2: Percolation thresholds for the 3D 0(2) model at different values of the field. 
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Figure 7: Kertesz line and pseudocritical line for the 3D 0(2) model. 



Again, the straight line in the log-log plot clearly indicates the power law behaviour 
of the Kertesz line. The slope of the straight line is 0.601(3), in good agreement with the 
exponent of the pseudocritical line 1/ (35 = 0.5995(6), obtained by using the 0(2) thermal 
exponents values determined in Hi~8[l . 

To compare the two lines directly, we make use of the results of a recent numerical 
determination of the 0(2) pseudocritical line [l~9[j. The results are included in Fig. |7|. The 
two lines are parallel in the log-log plot, since the exponents of h agree, but they do not 
coincide. 

We have shown that the Kertesz lines of the 2D Ising model and the 3D 0(2) model are 
in the limit of small fields described by the same functions as the corresponding thermal 
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pseudocritical lines, i.e., a power law with exponent 1/(58. This feature, which is likely to 
be general, suggests a relationship between the geometrical percolation singularities which 
still occur for H^O and the remaining non-singular thermal properties of the system. In 
both models we have found that the Kertesz line does not coincide with the pseudocritical 
line even in the scaling region we have explored. The difference in the behaviour of the 
two lines could well be a consequence of using the conventional h = bond weight (fj) also 
for h ^ 0; introducing a dependence on the field h into the Coniglio-Klein factor might 
well make the two lines coincide. Simple expressions of the form 1 — exp[— 2/3(l + h)] would 
lead to a Kertesz line which goes to infinity for h — > oo, so that an eventual interpretation 
of the pseudocritical line as a curve of percolation points is not yet excluded. Further 
work on this is in progress. 
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